The anisotropic Heisenberg model in a transverse magnetic field 
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One dimensional spin-1/2 XX Z model in a transverse magnetic field is studied. It is shown that 
the field induces the gap in the spectrum of the model with easy-plain anisotropy. Using conformal 
invariance the field dependence of the gap at small fields is found. The ground state phase diagram 
is obtained. It contains four phases with different types of the long range order (LRO) and a 
disordered one. These phases are separated by critical lines, where the gap and the long range order 
vanish. Using scaling estimations and a mean-field approach as well as numerical calculations in 
the vicinity of all critical lines we found the critical exponents of the gap and the LRO. It is shown 
that transition line between the ordered and disordered phases belongs to the universality class of 
the transverse Ising model. 



I 

o 
o 



> 
o 
o\ 
in 

o 

(N 
O 



o 
o 



X 



75.10.Jm - Quantized spin models 

The effect of a magnetic field on an antiferromagnetic 
chain has been attracting much interest from theoretical 
and experimental points of view. In particular, the strong 
dependence of the properties of quasi-one-dimensional 
anisotropic antiferromagnets on the field orientation was 
experimentally observed [1] . So, it is interesting to study 
the dependence of properties of the one-dimensional an- 
tiferromagnet on the direction of the applied field. The 
simplest model of the one-dimensional anisotropic anti- 
ferromagnet is the spin-1/2 XXZ model. This model in 
an uniform longitudinal magnetic field (along the Z axis) 
was studied in a great details [2]. Since the longitudinal 
field commutes with the XXZ Hamiltonian the model 
can be exactly solved by the Bethe ansatz. This is not 
the case when the symmetry-breaking transverse mag- 
netic field is applied and the exact integrability is lost. 
Because of a mathematical complexity of this model it 
has not been studied so much. From this point of view it 
is of a particular interest to study the ground state prop- 
erties of the ID XXZ model in the transverse magnetic 
field. The Hamiltonian of this model reads 
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(1) 



with periodic boundary conditions and even TV. 

The spectrum of the XXZ model for —1 < A < 1 is 
gapless. In the longitudinal field the spectrum remains 
gapless if the field does not exceed a saturated value 1-l-A. 
On the contrary, when the transverse magnetic field is 
applied a gap in the excitation spectrum seems to open 
up. It is supposed [3] that this effect can explain the pe- 
culiarity of low temperature specific heat in Yb4 Asa [1] . 
The magnetic properties of this compound is described 
by XXZ Hamiltonian with A w 0.98 and it was shown 
[3] that the magnetic field in easy plain induces a gap 
in the spectrum resulting in a dramatic decrease of the 
linear term in the specific heat. 



First of all, what do we know about the model (1)? 

The first part of the Hamiltonian is well-known XXZ 
model with the exact solution given by Bethe ansatz. In 
the Ising- like region A > 1 the ground state of XXZ 
model has a Neel long-range order (LRO) along Z axis 
and there is a gap in excitation spectrum. In the region 
— 1 < A < 1 the system is in the so-called spin- liquid 
phase with a power-law decay of correlations and a linear 
spectrum. And, finally, for A < — 1 classical ferromag- 
netic state is a ground state of XXZ model with a gap 
over ferromagnetic state. 

In the transverse magnetic field the total is not a 
good quantum number and the model is essentially com- 
plicated, because the transverse field breaks rotational 
symmetry in X — Y plane and destroys the integrabil- 
ity of the XXZ model, except some special points. In 
particular, the exact diagonalization study of this model 
is difficult for finite systems because of non-monotonic 
behavior of energy levels. 

The first special case of the model (1) is the limit 
A ±oo. In this case the model (1) reduces to the 
ID Ising model in a transverse field (ITF), which can be 
solved exactly by transformation to the system of non- 
interacting fermions. In both limits the system has a 
phase transition point he — | A|/2, where the gap is closed 
and the LRO in Z direction vanishes. 

It is suggested [5] that the phase transition of the ITF 
type takes place for any A > at some critical value h = 
hc{A). One can expect also that such a transition exists 
for any A and the transition line connects two limiting 
points he — |A|/2, A ±oo. 

Similar to these limiting cases, for any |A| > 1 and 
h < /ic(A) the system has the LRO in Z direction (Neel 
for A > 1 and ferromagnetic for A < —1). But for |A| < 
1 and h < hc{A) the ground state changes and instead of 
the LRO in Z direction a staggered magnetization along 
Y axis appears at /i < hc{A). 

This assumption is confirmed on the 'classical' line 
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FIG. 1. Phase diagram of the model (1). The thick sohd 
hnes denote the critical lines, thin sohd hne is the 'classical' 
line, and dashed line denotes the line /ii(A) (see below). 



/ici = ^2(1 + A) (/ici < ft.c(A)), where the quantum fluc- 
tuations of WZ model are compensated by the trans- 
verse field and the exact ground state of (1) at ft. = /ici 
is a classical one [4]. The excited states on the classical 
line are generally unknown, though it is assumed that 
the spectrum is gapped. 

The second case, where the model (1) remains inte- 
grable, is the isotropic AF case A = 1. In this case 
the direction of the magnetic field is not important and 
the ground state of the system remains spin-liquid one 
up to the point h = 2, where the phase transition of the 
Pokrovsky-Talapov type takes place and the ground state 
becomes completely ordered ferromagnetic state. 

And the last special case is A = — 1. In this case the 
model (1) reduces to the isotropic ferromagnetic model in 
a staggered magnetic field. This model is non-integrable, 
but as was shown [13], the system remains gapless up to 
some critical value h — ho, where the phase transition of 
the Kosterlitz-Thouless type takes place. 

Summarizing all above, we expect that the phase dia- 
gram of the model (1) (on (A, h) plane) has a form as 
shown on Fig.l. The phase diagram contains four regions 
corresponding to different phases and separated by tran- 
sition lines. Each phase is characterized by its own type 
of the LRO: the Neel order along Z axis in the region 

(1) ; the ferromagnetic order along Z axis in the region 

(2) ; the Neel order along Y axis in the region (3); and 
in the region (4) there is no LRO except magnetization 
along the field direction X (which, certainly, exists in 
all above regions). Hereafter under LRO we mean the 
corresponding to given region type of LRO. 

In this paper we investigate the behavior of the gap and 
the LRO near the transition (critical) lines. The first sec- 
tion is devoted to the classical line, where we review the 
exact ground state and construct three exact excitations. 
In the section 2 we study the transition line hc{A) with 
use of the mean-field approach and exact diagonalization 
of finite systems. In the section 3 we find the critical 
exponents in the vicinity of the line h — 0. The proper- 



ties of the model near the critical lines A = ±1 and in 
particular in the vicinity of the points (A = ±1, h — 0) 
are studied in sections 4 and 5. 



I. THE CLASSICAL LINE 



First we consider the classical line hc\ — ■\/2(l -|- A) 
(A > —1), because further we shall often refer to this 
line. This line is remarkable in a sense that the ground 
state on it is identical to the classical one and quantum 
fluctuation are missing. It was shown in [4] that the 
ground state of (1) on this line is two-fold degenerated 
and the ground state wave functions with momentum 
k = and k = n are 

and $1(2) are direct product of single-site functions: 
|$i) = \a1a2a3a4 . . .) 

|$2) = \a1a2Oi3a4 . . .) 

where \ai) is the state of i-th spin lying in the XY plane 
for |A| < 1 (or in XZ plane for A > 1) and forming an 
angle (p with X axis. These states can be written as: 

K) = (e''^5+ - 1) li) , |A|<1 
\a,) = (e^5+ - 1) li) , A > 1 

with cos 93 — hc\/2 for |A| < 1 and cosh ip — hc\/2 for 
A > 1. 

The state \ai) is obtained by rotation of i-th spin by tt 
about the axis of the magnetic field X. 

|a,) = e'^^? |a,) 

The ground state has a two-sublattice structure and 
is characterized by the presence of the LRO in the Y 
(|A| < 1) or in the Z {A > 1) directions. In particular, 
for |A| < 1 the staggered magnetization (5*1^) is 

(SI) = ^^/^^ 

However, the excited states of (1) on the classical line 
are nontrivial in general. Nevertheless, some of them 
can be found exactly. For this aim it is convenient to 
introduce the operator overturning the i-th spin: 

i?. = e'"^-, |A|<1 
R, = e'''^", A>1 

so that the states of the 'overturned' i-th spin \f]i) = 
Ri lofi), \l3i) = Ri \ai) are orthogonal to \ai), \ai): 

(a, I A) = (a,| A) =0 
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Then, the exact excited state are written as: 



(2) 



'tn 

= ^(-iri?„i?„+il$i(2)) 

n 

-^(-iri?„i?„+ii?„l$i(2)) 



So, each of three exact excitations is also two-fold de- 
generated. Actually, this degeneracy is a consequence of 
the Z2 symmetry describing the rotation of all spins by 
TT about the axis of the magnetic field X. 

To show that these states are really exact ones it is 
convenient to rotate the coordinate system, so that one 
of the groimd states, for example $1, will be all spins 
pointing down. For the case |A| < 1 this transformation 
is the rotation of the spins on even (odd) sites by an 
angle ip {—(fi) around the Z axis followed by the rotation 
by 7r/2 about the Y axis: 



= COS ip + {-ly al sin ifi 
SI = {-iTcrlsmip - alcos^p 



QZ —X 



(2) 



For A > 1 the transformation of the spin operators is 
defined by the relation 



SI = al cos + (- sin ^ 



Sn — '^n 



5^ = -(-!)"< sin + < cos ^ 



(3) 



Then the Hamiltonian (1) on the classical line becomes 
ifi = A^a„-a„+i + (l + A)^< 



+^ciV 1 - -f + <-i + 1) (4) 

for A < 1 and 

H2 = Y,^n- <Tn+l - (A - + 



+A'i-4E(-i)"<«+i + <-i + 1) (5) 

for A > 1. 

The ground state of both Hamiltonians as well as of 
(1) is two-fold degenerated. One of the ground state is 
obviously all spins cr„ pointing down $1 = |0) = ...). 
The energy of this state is: 



N ,,A 

Eo = N— 

" 2 4 



(6) 



The second ground state $2 in this representation has 
a more complicated form: 

$2 =n(cos^ + (-l)V+sin(^)|0) 



Now, it is easy to see that the following three excited 

states are exact ones: 



(1) 



(2) 



(3) 



^a+|0), E,-Eo = l + A (7) 
^(-l)V+<7++i|0), £;2-£;o = 2 + A 

V+<7++ia+ \0),Es-Eo = 3 + 2A 



One can check that last terms in (4) and (5) anni- 
hilate these three functions and, hence, they are exact 

excited states of (1) for any even A^. Similarly to the 
ground state the excited states (7) are degenerated with 
the states |V'2)- These states j^/ij) can be represented in 
the same form (7), but in the coordinate system, where 
the function ^2 is all spins pointing down. 

The states V'i(2)^ especially interesting because 

they define the gap of the model (1) on the classical line 
at small value of hc\. Our numerical calculations of fi- 
nite systems show that at hd ^ (A ^ —1) the lowest 
branch of excitations has a minimum at k = and corre- 
sponding excitation energy is (1 -|- A) (of course, due to 
the Z2 symmetry there is another branch with the min- 
imum at k = TT and the same minimal energy, but we 
consider one branch only). Excitation energy at fc = tt 
obtained by the extrapolation of numerical calculations 
at iV — > 00 is 2(1 + A). When hc\ increases excitation 
energies at A; = and k = tt are drawn together and at 
some hc\ they are equal to each other. Our numerical 
results give hd ~ 0.76 (A ~ —0.79). So, the gap on the 
classical line is (1 -|- A) for -1 < A < -0.79. 



II. THE TRANSITION LINE H = He {A) 

The existence of the transition line /ic(A) passing 
through the whole phase diagram is quite natural, be- 
cause at some value of the magnetic field all types of the 
LRO except the LRO along the field must vanish. The 
transition line connects two obvious limits A — > ±00, 
when the model (1) reduces to the ITF model. The line 
passes through the exactly solvable point (A = 1, /i = 2) 
and the point (A = —1, h = ho) studied in [13]. We 
suppose that the whole line hc{A) is of ITF type with 
algebraically decaying of correlations with corresponding 
critical exponents [15]. 

The transition line can be also observed from the nu- 
merical calculations of finite systems. As an example, the 
dependencies of excitation energies of three lowest levels 
on h and for = 10 — 18 are shown on Fig. 2. From this 
figure one can see that two lowest states cross each other 
N/2 times and the last crossing occurs on the classical 
point hci = V2. These two states form two-fold degen- 
erated ground state in the thermodynamic limit. They 
have different momenta k = and fc = tt and different 
quantum number describing the Z2 symmetry, which re- 
mains in the system after applying the field. As for the 
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FIG. 2. The dependence of difference of tlie energy of two 
lowest levels and the ground state energy on magnetic field h 
for finite chains with iV = 10, . . . 18. 



first excitation above the degenerated ground state, on 
Fig. 2 we see also numerous level crossings. These level 
crossings lead to the incommensurate effects which man- 
ifest itself in the oscillatory behavior of the spin correla- 
tion functions. The correlation functions at n l^- 1 have 
a form: 



(5f 5^) - = /(n)e- 



(8) 



where (S"") {a = x, y, z) is the corresponding magnetiza- 
tion (the LRO), and f{n) is the oscillatory function of n 
with the oscillating period depending on h and A. All 
crossings disappear at h > /ici(A) and in this region of 
the phase diagram the correlation functions do not con- 
tain oscillatory terms. 

The energy of the first excitation near h^i goes down 
rapidly, and after extrapolation we found that for the 
case A = at the magnetic field he — 1.456(6) > hc\ the 
gap vanishes. Inside the region hd < h < he the ground 
state remains two-fold degenerated, though there is no 
level crossings. At h > he the mass gap appears again 
and for a large field the gap is proportional to h. 

In order to determine the transition line /ic(A) and 
to study the model in the vicinity of hc{A), we use the 
Fermi-representation of (1). This representation gives 
the exact solution in the limits A — > ±oo and, besides, it 
yields the exact ground state on the classical line. 

At first it is convenient to perform in (1) a rotation of 
the spins around the Y axis by 7r/2, so that the magnetic 
field is along the Z axis: 

if = ^(A5^5^+i + 5^5^+i + 5^5^+i) + /i^5^ (9) 

After Jordan- Wigner transformation to Fermi opera- 
tors a+, a„ 



Si 



1 
2 



a' a 



n 



(10) 



the Hamiltonian (9) takes the form 



hN N 
1- A 



{h-1 



1 + A 



cos k)a^ak 



4 ^smk{a+atk+a-kak) 



(11) 



Treating the Hamiltonian Hf in the mean-field approx- 
imation we find the ground state energy Eq and the one- 
particle excitation spectrum e{k): 



i^o/iV = (.-l)(7t-i)+i-(l-f) 



72 



+§73 +7? -71+73 



(12) 
(13) 



e{k) = ^/a^{k) + b^{k) 
where g = 1 — A and 

a{k) = (/i - 1) - (^1 - 1^ cos k + 2ji- 272 cos k 

6(fc) = (|+273)sinfc (14) 

Quantities 71 , 72 and 73 are the ground state averages, 
which are determined by the self-consistent equations: 



71 



72 = (a+a„+i) 



fc>0 



a{k) 



73 = (aia+^i) = -Yl 



e{k) 

E^T^ cos A; 
fc>o ^ 

h{k) 



fe>0 



2e{k) 



sin k 



(15) 



Magnetization S = (S^) of the model (1) is equal to 



5=--7t 



(16) 



The numerical solution of Eqs.(15) shows that the 
function e(fc) has a minimum at fcmin, which is changed 
from 7r/2 at ft, = to zero at h = /ii(A) and fc,„in ~ 
for h > hi{A). The gap in the spectrum e(fc) vanishes at 
he{A) {he > hi) and for /i > /ii is m = |/i — /id- The de- 
pendencies of hi{A) and hc(A) are shown on Fig.l. We 
note that the Hamiltonian H{ differs from the domain- 
wall fermionic Hamiltonian which is mapped of (1) in [5]. 
The transition line obtained in [5] is a linear function of 
A in contrast to /ic(A) on Fig.l. 

It is interesting to note that the mean field approxi- 
mation gives the exact ground state on the classical line 
hci — \/2(l -I- A). On this line the solution of Eqs.(15) 
has very simple form: 



1 hc\ 

71 = 2 - 



72 



'73 



71 



1 



hi 

72 = 73 = — 



16 
-4 



\A\ < 1 



A > 1 (17) 



4 



and the energy Eo/N = — i — ^. 

The gap on the classical line in the mean-field approx- 
imation is 



m 



|A| <1 



m=42^(/ici-2)2, A>1 



2h 



(18) 



cl 



We compared (18) with the results of the extrapolation 
of finite systems on the classical line. The coincidence is 
fairly good for A > 0.5. Eq.(18) gives rather satisfactory 
estimation for the gap up to A ~ —0.5. For example, 
at A = (ftci = \/2) from Eq.(18) m = 0.086, while 
extrapolated gap is w ~ 0.076(4). 

The smaller the fermion density the better works the 
mean-field approximation. It becomes worse when the 
magnetization S ^ 0. This is the reason of incorrect 
behavior of the gap at /id ^ (A ^ —1). According to 
(18) m = 1, while it vanishes in this limit as m = 1 -|- A 
(7). 

The Hamiltonian Hf in the mean-field approximation 
has a form similar to the well-known bilinear Fermi- 
Hamiltonian describing the anisotropic XY model or the 
ITF model. Using results of [15] the following facts re- 
lated to the considered model can be established: 

1. There is a staggered magnetization {S^) along the 
Y axis for |A| < 1 or (5^) along the Z axis for |A| > 1 
and they vanish at h ^ he as {he — hY^^. 

2. The magnetization S has a logarithmic singularity 
at h ^ he- 

3. The spin correlation function decay exponentially 
(excluding the transition line) at n — > oo: 



G«(n) = (5f 5«) - = /(n)e- 



(19) 



The function /(n) has an oscillatory behavior at < /i < 
hci, while at h > hd it is monotonic. At h = hd, f{n) = 



and f{n) 



at {hci-h) < 1. Thus 



the classical line determines the boundary on the phase 
diagram where the spin correlation functions show the 
incommensurate behavior. 

On the transition line h = hc{A) the spin correlation 
functions have power-law decay: 



G^(n) ~ l/r^^ 
for |A| < 1 and 



h^ = 


2 32' 






hi = 


hr — —, 

" 16' 






m = 


\h — hc\, h > hi 






m = 


2s/2\ 32' 


h<hi 


(20) 



The magnetization S is 
S=\- —^/hc-h-^, g^^hc-h 



9_ 
4tt 



2(fec - h) 
-Kg 



In 



9 



hr — h 



^VK-h (21) 



The susceptibility x(/i) = ^ is 



X{h) = - In f ^ 
ng \hc- 

Xih) ^ ^ 



, 5 > - h 

g «; ^/hc-h (22) 



\/2tt \Jhc — /i' 

As follows from (22) there is a crossover from square 
root to logarithmic divergence of x when the parameter 

2 

■j^rih changed from to oo. 



III. THE LINE = 0, |A| < 1 

A. Scaling estimations 

The XXZ model is integrable and its low-energy prop- 
erties are described by a free massless boson field theory 
with the Hamiltonian 

Ho = lJdx [tf + {d,^f] (23) 

where n(x) is the momentum conjugate to the boson field 
$(a;), which can be separated into left and right moving 
terms $ = $l -|- The dual field $ is defined as a 
difference $ = $i — The spin-density operators are 
represented as 

5^^^a,$ + const.(-l)"cos| 



Gy{n) ~ l/r^3/^ G'{n) ~ l/n'^\ |Af^^ cos hnm 



C(-l)" + const, cos — 



(24) 



for |A| > 1. 

These results show that the transition at h = hc{A) 
belongs to the universality class of the ITF model. 

In the vicinity of the point /i = 2, A = 1 the fermion 
density is small (S* ~ i) and the mean-field approxima- 
tion of the four fermion term gives the accuracy, at least, 
up to g^ or (2 — h^. For this case we give corresponding 
expressions (at 5 <C 1): 



with constant C found in [7] . The compactification radius 
R is known from the exact solution 



2nR^ = 6 



arccos(A) 



The non-oscillating part of the operator in Eq.(24) 
has scaling dimension d = 9/2 + 1/29 and confer mal spin 
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S = 1. The non-zero conformal spin of the perturbation 
operator can cause the incommensurabihty in the sys- 
tem [8], which is in accord with Eq.(19). As was shown 
in [9] , the common formula for the mass gap 



1 



2-d 



1/0 



(25) 



is not apphcable in the whole region |A| < 1. Due to 
non-zero conformal spin of the non-oscillating part of 
the operator it is necessary to consider higher-order 
effects in h. The analysis shows [9] that the original per- 
turbation with nonzero conformal spin generates another 
perturbation with zero conformal spin 



V = h'^ cos (47ri?$) 



(26) 



This perturbation gives the critical exponent for the 
mass gap 



7 



1-61 



(27) 



Comparing Eq.(25) and (27) we see that the pertur- 
bation (26) becomes more relevant in the region A < 
cos[ttV2] « -0.266. 

It turns out that the oscillating part of the operator 
gives another, more relevant index for the gap at A < 0. 
Let us reproduce usual conformal line of arguments for 
this oscillating part. The perturbed action for the model 



IS 



S = Si 



dtdxS''{x,t) 



(28) 



where So is the Gaussian action of XXZ model. The 

time-dependent correlation functions of the XXZ chain 
show the power-law decay at | A| < 1 and have the asymp- 
totic form [6] 



{S^{x,r)S^O,Q)) 



(-l)"Ai 

(29) 



with Ai, A2 are known constants [7] and r = it is imagi- 
nary time. Therefore, we can estimate the large-distance 
contribution to the action of the oscillating part of the 
operator S^{x,t) as 



h j dTdxS''{x,T) 



hjdrY, 



{-If 



n=2l 



9x 



h / drdx 



(a;2+t;2r2)3+i 



The relevant field S'^(a;,r) gives rise to a finite corre- 
lation length ^. This correlation length is such that the 
contribution of the field 5"^(a;, r) to the action is of order 
of 1. That is 
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7 on A. The smallest exponent gives the most relevant type 
of perturbation and defines the index for the mass gap. 



/ 

Jo Jo 



dx 



Ox 



(3.2 + ^2^2)9/4-M 



which gives for the mass gap 



6/2 



(30) 



In fact, the oscillating factor (—1)" in the correlator 

in some sense eliminates one singular integration over x, 
and into common conformal formula m ~ h'^-'^ , where D 
is the dimension of space and d is the scaling dimension 
of perturbation operator, one should put D = 1 instead 
of conventional D = 2. 

The comparison of the expressions Eqs.(25), (27) and 
(30) shows that for < A < 1 the leading term is given 
by Eq.(25), while for -1< A < by Eq.(30). Thus, one 
has: 



< A < 1 
-1< A < 



(31) 



The functions z^(A), /i(A), 7(A) are shown on Fig.3. 

In this respect the model (1) is different from the XXZ 
model in the staggered transverse field for which m ~ 

f^2/i4-0) for all < 1 [10]. 

The staggered magnetization (LRO) along Y axis be- 
haves as: 



(32) 



Hence, the LRO has also two different critical expo- 
nents: 



{\sy\) - h'^-o-i/o < A < 1 
{\sy\) ~ - 1< A < 



(33) 
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B. Perturbation series 

The critical exponents v and \x can be also derived 
from the analysis of infrared divergences of the perturba- 
tion theory in h. Obviously only even orders in h gives 
contribution. Let us estimate the large-distance behav- 
ior of the following operator, defining the order of the 
perturbation series: 



V = 



1 



Eq — Hq Eq — Hi 



-V 



(34) 



with V = hY^Sf and Hq is the Hamiltonian of XX Z 
model. The perturbation series for the ground state en- 
ergy has a form: 



1 



Eq — Ho 



V{l + U + U'^ + .. 



(35) 



Hereafter we consider large but finite system of the 
length N. Wo shall care about indices at N and h only, 
omitting all factors. 

At first we consider non-oscillating part of the corre- 
lator (29). Taking into account only low- lying excitation 
of linear spectrum of the XXZ model, giving the most 
divergent part, and estimating large-distance behavior of 
non-oscillating part of the correlator (29), we arrive at 



U 



{1/Nf 



lSfO+i/0 



A-e-i/e 



(36) 



Now we see, that if 4 — 6* — 1/0 > 0, then each next 
order in perturbation theory (35) diverges more and more 
strongly. To absorb these infrared divergences one has to 
introduce the scaling parameter y ~ Nh'^ and guess that 
the series (l + U + U"^ + . . .) in (35) forms some function 
of the scaling parameter 2/. In our case = 2/(4—0—1/0) 
(see Eq.(25)) and U ~ y^/^^ 

The leading divergence of the second order of the 
ground state energy can be found in a similar way as 
in (36): 



-V ~ /i27v3-e-i/« 



(37) 



Eq — Ho 

Combining the formulae (36) and (37), we can rewrite 

SE ~ Nh^'-fiy) 

with some unknown function f{y) having the expansion 
at small y: 

In the thermodynamic limit N ^ oo the scaling pa- 
rameter y = Nhy also tends to infinity y ^ oo. Since 
the energy is proportional to N, the function f{y) has a 



finite limit /(oo) = a. Thus, in the thermodynamic limit 
for the correction to the ground state energy one has: 



5E ~ aNK 



2v 



(38) 



For the first excitation state the divergences of the or- 
ders of the perturbation theory have the same form as in 
(36) and (37). So, for the gap one finds the same scaling 
parameter y = Nh" and 

In the thermodynamic limit the mass gap is of order 
of unity (in terms of A^), and, therefore, the function 
g{y) ~ Xjy at y — > oo. Thus, finally we arrive at the 
Eq.(25). 

Now we consider more subtle, oscillating part of the 
correlator (29). For the oscillating part at large distances 
one can write: 



(-ir 



1 



1 



It turns out that the oscillating part of the pertur- 
bation operator V connects the low-lying gapless states 
with the finite energy states. That is, each second level 
in all orders of the perturbation series has a finite gap 
from the ground state. Therefore, for the operator U one 
has 



^2 k^i^f) 



(l/AT) 



Since Q is always less than 2 the divergences grow with 
the order of the perturbation theory. To eliminate these 
divergences we introduce the scaling parameter y = Nh^, 
with ^ defined in Eq.(30), so that U ^ y^^^. 

The second order for the ground state energy in this 
case is 

/qxcx\ 

and the total correction to the ground state energy has 
a form 

5E ~ Nh^'^^f{y) 

with some unknown function f{y), having the finite limit 
/(oo) = b. 

Thus, in the thermodynamic limit the ground state 
energy behaves as 

SE ~ bNh^i^ 
The mass gap is found similarly 

TO ~ Nh^fgiy) 

with the function g{y) ~ 1/y aX y ^ oo. Thus, in the 
thermodynamic limit we reproduce Eq.(30). 
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We note that wc have estimated only long wave-length 
divergent part of the perturbation theory. Besides, there 
is a regular part of the perturbation theory, which gives 
the leading term ~ /i^. So, combining all the above facts 
we arrive at 

SE 



N 2 



(39) 



As it can be seen from Eq.(39), 5E consists of a regular 
term and two singular terms. Since v > 1 and /i > 
1, the susceptibility x is finite for any A in contrast to 
the model with staggered transverse field [10], where the 
singular term is with rj = A/ [A — 9) < 2. 

From Eq.(25) and (30) it follows that 1 at A ^ 1 
and /U^latA^— 1. Hence, in both limits one of the 
singular terms becomes proportional to , and, there- 
fore, gives contribution to the susceptibility. It implies 
that in the symmetric points A = ±1 the susceptibility 
has a jump. 



IV. THE LINE A = 1 

In the vicinity of the line A = 1 it is convenient to 
rewrite the Hamiltonian (1) in the form 

H = Hq + V 

Ho = ^ S„ • S„+i + h ^ 



(40) 



where the parameter (7=1 — A<Clis small. On the 
isotropic line A = 1 the model (1) is exactly solvable by 
Bethe ansatz. The properties of the system remain crit- 
ical up to the transition point he = 2, where the ground 
state becomes ferromagnetic. Therefore, for /i < 2 and 
small perturbation V we can use conformal estimations. 
The asymptotic of the correlation function on this line 

is 



(-1)" 



(41) 



where a{h) is known function obtained from Bethe ansatz 
[11] and having the following limits: 



a{h) 



a{h) - -, 



1 



21n(l//i)' 

h^2 



(42) 



So, the scaling dimension of operator 5*^ is = 
a{h)/2, and the scaling dimension of operator SfSf^^ 
is four times greater d^z = Ad^ = 2a{h). Since a{h) < 1, 
then the perturbation V is relevant and leads to the mass 
gap and the staggered magnetization 



m 



{\S'\) 



\9\ 
\9\ 



l/(2-ci..) ^ |^|l/(2-2a) 



A < 1 
A > 1 



a/(4— 4a) 
|a/(4— 4a) 



(43) 



From the general expressions for the mass gap (43) , in 
the limit h ^ 2 we obtain that m ^ g, which is in accord 
with the result of the mean field approximation (20). 

The LRO in the vicinity of the point A = 1, h = 2 
vanishes on both lines: at A = 1 from (43) as g^^^; and 

a,t h = he as \hc — h\^^^. Besides, one has the exact 
expression for LRO on the classical line 



2V2 



(44) 



Combining all these facts we arrive at the following for- 
mula: 



{\Sy\)=2-y^g'^^\he-h\'/' 



(45) 



The behavior of the system near the point A = 1, 
/i = is more subtle. As it follows from Eq.(31), for very 
small h the mass gap is m ~ /i, while on the other hand 
from Eq.(43) one obtains another scaling m --^ g^'^('^/'^\ 
So, there are two regions near this point with different 
behavior of mass gap. The boundary between these two 
regions can be found from the following consideration. 
Let us rewrite the perturbation in the Hamiltonian (40) 
in the form: 

V = Vi + V2 

^2=2 ~ ^n^n+l) 

The part of the Hamiltonian Hq + Vi corresponds to 
the XX Z model in the longitudinal magnetic field, which 

is gapless when the magnetic field h > exp (— 7r^/2,y^). 
Therefore, in the region of very small magnetic field h < 
exp (—• 7r^/2y^) the perturbation Vi is relevant, leading 
to the mass gap m ^ h. The two-cutoff scaling procedure 
[8] [9] gives rise to the mass gap m ~ /lexp (— 7r^/2,yg) 
for h > exp (— 7r^/2y^) . And, finally, when g is much 
less than the scaling dimension of opcTator V2 defines 
the exponent for the gap (43). Summarizing all above, 
the mass gap in the vicinity of the isotropic point A = 1, 
/i = is: 



m h, 



g 



-Inh 



Inh < — ^ 

—J- » In/i » -4= 

ln/i»— i (46) 



V. THE LINE A = -1 

In this section we consider the model (1) in the vicinity 
of the line A = —1, where (1 + A) = 5 <C 1 is a small 
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parameter. It is convenient to rotate spins on each odd 
site by tt around Z axis, so that the model (1) becomes: 

F = - ^ S„ • S„+i + 8 ^ 5^5^+1 - h Y,{-\TSl 

(47) 

At (5 = and ft = the ground state of (47) is fer- 
romagnetic state (degenerated with respect to total S^) 
with zero momentum. The states which can be reached 
from the ground state by means of the transition operator 
^(— 1)"'S'^ are the states with g = tt and finite gap over 
the ground state. When ^ ^ 1 the transition operator 
connects the low-energy states and the states with ener- 
gies £s ~ 2. The second order correction to low-energy 
states is 



8Ef^ = ft^ E 



(;|(-i)"5-|g) (s| (-1)-^- 10 

El - Es 



(48) 



where \l) is the low-energy state and \s) is the state with 
high energy ~ 2. So, for 5 <C 1 Eq.(48) can be rewritten 
as 



(49) 



nm 



The spin correlation function (/| S^S^ \ l) is slow vary- 
ing function of \m — n\ &t 5 ^ 1. Therefore, 

{l\ {-IT-'SISI, \l) ('I ^n^n+l 10 (50) 



Thus, as follows from Eqs. (48-50), the low- lying states 
of (47) at \5\ < 1 and /i «; 1 are described by the XYZ 
Hamiltonian 

(51) 

The coincidence of the low-energy spectra of (47) and 
(51) in the vicinity of ferromagnetic point A = — 1, /i = 
has been checked numerically for finite systems. The 
spectrum of low-lying excitations of the s — ^ XY Z 
model (51) as well as initial model (1) near the ferro- 
magnetic point A = — 1, /i = can be described asymp- 
totically exactly by the spin-wave theory, which gives 



A > -1 



TO = /i^(l + A)/2, 
m= V(l + A)(l + A-}-ft2/2), 



(52) 



A < -1 



It can be checked that Eq.(52) yields the exact gap of 
the XYZ model [14] at \5\, /i < 1. The validity of the 
spin-wave approximation is quite natural because in the 



vicinity of the ferromagnetic point A = — 1, ft = the 
number of magnons forming the ground state is small. 

We note also that the gap (52) for A > —1 is in accord 
with the conformal theory result (31) and gives us the 
preexpo nential fac tor for the gap. On the classical line 
ftci = \/2(l-|- A) Eq.(52) yields the gap m = 1 -|- A 

which confirms that the function tj)^^ in Eq.(7) gives the 
exact gap. 

The similar mapping of the model (1) with arbitrary 
spin s to XYZ model can be done for A ~ —1, ft -C 1. 
Taking into account that Eg = 4s corresponding XYZ 
Hamiltonian is 



H 



[(1 



+ 1 



cy qV 
'-'n'-'n+l 



AS'^S'^+i] 



■^E(^n)^ (53) 



where S" are spin-s operators. 

The leading term of the gap of the model (1) with 
arbitrary spin s in the vicinity of the point A = — 1, 
ft = is given exactly by the spin-wave theory: 



A > -1 



m = ftA/(l A)/2, 
TO = 2sV(l + A)(l-FA-hftV8s2), 



(54) 



A < -1 



On the classical line ftci = 2s\/2{l + A) Eq.(54) gives 
the correct result to = 2s5. 

Strictly on the line A = — 1 the model (1) reduces to 
the isotropic ferromagnet in the staggered magnetic field. 
This model is non-integrablc. but it was suggested [13], 
that the system is governed by a c = 1 conformal field 
theory up to some critical value ft = fto, where the phase 
transition of the Kosterlitz-Thouless type takes place. 

At ft < 1 where the mapping of (47) to the XYZ 
model (51) is valid, the line A = — 1 is described by the 
XX Z model and the correlation functions have a power- 
law decay: 



/ Qz \ — / qy qy \ V ' 



I QX QX \ 



(-1)" 

„B{h) 



(55) 



We believe that the relation between indices of X and 
y, Z correlators on the line A = — 1 has the form (55) for 
< ft < ftg. So, the scahng dimensions of operators Sf 
and Sf , on this line are dx = /3/2 and dy = dz = 1/2/3. 

On the fine A = — 1 the model (1) is gapless for h < ho. 
It means that the magnetic field term is irrelevant at 
ft < fto > 4) and becomes marginal at ft = fto, 

where dx = 2 and /3(fto) = 4. So, at the point ft = fto 
the transition is of the Kosterlitz-Thouless type, and for 
ft > fto the mass gap is exponentially small. 

In the vicinity of the line A = — 1 the term S J2 ^nSn+i 
in (47) can be considered as a perturbation and the scal- 
ing dimension of the perturbation operator S^S^_^_l is 
dzz = Mz = 2//3(ft). Since /3(ft) > 4 at ft < fto, the 
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perturbation is relevant and leads to the mass gap and 
the LRO: 



m ' 

{\sy\) 
i\s'\) ■ 



,2/3 



(5 > 
(5 < 



(56) 



In particular, m ^ \5\''^ and (l^^'l) ~ 1^1'/^^ at /i ^ 
ho. 

The function /3(/i) is generally unknown, except the 

case /i <C 1, where the mapping to XX Z model is valid, 

and (3{h) = 1 — iarccos(-^2 — 1) — T^/h. However, 

since the model (1) at A = —1 and h < ho is conformal 
invariant, we can use a finite-size scaling analysis to de- 
termine the exponent (3{h) and the vahie of ho- Accord- 
ing to the standard scaling approach [16] P{h) = 
where v is the sound velocity and ^ is the difference be- 
tween the two lowest energies of the system. We carried 
out calculations of for finite systems. The extrapo- 
lated function f3{h) agrees well with the dependence w/h 
at /i <C 1 and /3 = 4 at /iq — 0.52. This estimation is close 
to our direct numerical estimations ho — 0.549. On the 
other hand, the mean-field approach gives rather crude 
value ho = hc{—l) = 0.69. 



VI. CONCLUSION 

In summary, we have studied the effect of the 
symmetry-breaking transverse magnetic field on the s = 
1/2 XX Z chain. On the contrary to the longitudinal 
field the transverse field generates the staggered magne- 
tization in Y direction and the gap in the spectrum of 
the model with the easy-plain anisotropy. Using the con- 
formal invariance we have found the critical exponents 
of the field dependence of the gap and the LRO. It was 
shown that the spectrum of the model is gapped on the 
whole {h, A) plain except several critical lines, where the 
gap and the LRO vanish. The behavior of the gap and 
the LRO in the vicinity of the critical lines A = ±1 is 
considered on the base of the conformal field theory. We 
note that in the vicinity of the points (A = 1, /i = 0) 
and (A = 1, /i = 2) there is crossover between different 
regimes of the system behavior. It is shown that near the 
point (A = — 1, h = 0) the initial model can be mapped 
to the effective exactly solvable ID XYZ model and has 
the spin- wave spectrum. The transition line hc{A) be- 
tween the ordered phases and the disordered one is stud- 
ied in the mean-field approximation. This study shows 
that this transition is similar to that in the Ising model in 
the transverse field. However, the behavior of the model 
on the transition line near the Kosterlitz-Thouless point 
(A = — l,/i = /io)is not so clear. The mean-field approx- 
imation becomes worse at A — > — 1 and more sophisti- 
cated theory is needed. 
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